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for (y=0; y<M+3; ++ 

for <x=0; x<N+5; [ 

if ((x-3)>*l && (x-5)<=N-2 && (y-2)>=l && <y-3)<=M-2) { 
if ((x-5)>*l && (y-3)>=l) { 
if (out_compute ==255) { 

if (comp_edge_pixels [ ( (x-4)%3) *3+ (y-2) %3J <comp_edge_middle) out_compute=0 

if (comp_edge_pixels [ ( (x-4) %3) *3+ (y-4) %3] <comp_edge_middle) out_compute=0 

if (comp_edge_pixels I ( (x- 5) %3) *3+ (y-4) %3] <comp_edge_middle) out_compute=0 

) } 

if ((x-3)<=N-2 && (y-2)<=(M-2)) { 
maxdif f_compute = 

maxl3 (abs (gauss_xyjpixels [ ( (x-2) %3) *3+ (y-1) %3) 

- gauss_xy_middle) , maxdif f_compute) ; 

maxdif f_compute = 

maxl3 (abs (gauss_xyjpixels ( ( (x-2) %3) *3+ (y-3) %33 

- gauss_xy_middle) , maxdif f_compute) ; 

maxdif f_compute « 

maxl3 (abs (gauss_xy_pixels [ ( (x-3) %3) *3+ (y-3) %3) 

- gauss_xy_middle) , maxdif f_compute) ; 
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for (y*0; y<M+3; ++y) 
for (x«0; x<N+5; +4 

if (x>=4 && x<oN+3 && y>=3 && y<=M+l) { 
if((x-5)>»l && (y-3)>«l> I 
if (out compute 255) { 

csexmTn4mod3x3 » ((x-4)%3)*3; 
cseymin4mod3 * (y-4)%3; 

if (comp_edge_pixels [csexmin4mod3x3 + (y-2) %3] <comp_edge_middle) out_computesO ; 

if (comp_edge_j>ixels [csexmin4mod3x3 +cseymin4mod3] <comp_edge_middle) out_compute=0; 

if (comp_edge_pixels [ { (x-5)%3)*3 +cseymin4mod3] <comp_edge_middle) out_compute=*0; 

if ((x-3)<=N-2 && (y-2)<=(M-2)) { 
csexmin2mod3x3 = ( (x-2) %3) *3 ; 
cseymod3 = y%3; /* » (y-3)%3 */ 
maxdif f_compute • 

max 13 (abs (gauss_xy_pixels [csexmin2mod3x3 + (y-l)%3] 

- gauss_xy_middle) , maxdif f_compute) ; 

maxdif f_compute * 

max 13 (abs (gauss_xy_pixels tcsexmin2mod3x3 + cseymod3) 

- gauss_xy_middle) , maxdif f_compute) ; 

maxdif f_compute * 

maxl3 (abs (gaus_xy_pixels [ (x%3) *3 + cseymod3] 

- gauss_xy_middle) , maxdif f_compute) ; 
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distrib^ftrity: (x + 4)%3 = (x%3 ^V 3 > %3 

constant^xolding: = (x%3 %3 

constant unfolding: = (x%3 + 1%3)%3 

invert distributivity : = (x + 1)%3 

(a) 

modulo expansion: (x+2)%3 = 3 - x%3 - (x+l)%3 

(b) 
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for (y»0; y<M+3; ++y) 
cseymod3 « y%3; 
cseyminlmod3 ■ (y-1 
cseymin2mod3 <= (y-2)%3 
cseymin4mod3 « (y-4)%3; 
for (x=0; x<N+5; ++x) { 



# 



if (x>=4 && x<*N+3 && y>=3 && y<=M+l) { 
if((x-5)>=l &6 (y-3)>=l) | 
if (out compute »» 255) { 

csexmIn4mod3x3 * ( (x-4) %3) *3 ; 

if (comp_edge_pixels [csexmin4mod3x3 +cseymin2mod33 <comp_edge_middle) out_compute=0; 

if (comp_edge_pixels (csexmin4mod3x3 +cseymin4mod3] <comp_edge_middle) out_compute=0; 
if (comp_edge_pixels [ ( (x-5) %3) *3 +cseymin4mod3] <comp_edge_middle) out_compute~0; 

if ((x-3)<*N-2 && (y-2)<=(M-2) ) { 
csexmin2mod3x3 = ( (x-2) %3) *3 ; 
maxdiff compute = 

maxl3Tabs(gauss_xy_pixels [csexmin2mod3x3 + cseyminlmod3] 

- gauss_xy_middle) , maxdif f_compute) ; 

maxdiff compute = 

maxl37abs(gauss_xy_pixels [csexmin2mod3x3 + cseymod3] 

- gauss_xy_middle) , maxdif f_compute) ; 

maxdif f_compute ■> 

maxl3 (abs (gaus_xy_jpixels [ (x%3) *3 + cseymod3) 

- gauss_xy_middle) , maxdif f_compute) ; 
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cseymod3 » ^^^B ^^^B 

for (y=0; y<M+3; ++y) 

cseyminlmod3 » cseymod3; 

cseymod3 = y%3; 

cseymin2mod3 » 3-cseymod3-cseyminlmod3 ; 
for (x*0; x<N+5; ++x) { 

if (xi»4 && x<=N+3 && y>-3 &fc y<*M+l) { 
if((x-5)>«l t& (y-3)>=l) j 
if (out compute == 255) { 

csexmIn4mod3x3 - ({x-4)*3)*3; 

if (comp_edge_pixels [csexmin4mod3x3 +cseymin2mod3) <comp_edge_middle) ouc_compute«0 ,- 

if (comp_edge_pixels [csexmin4mod3x3 +cseyminlmod3] <comp_edge_middle) out_compute»0; 
if (comp_edgejpixels t ( (x-5) %3) *3 +cseyminlmod3 J <comp_edge_middle) out_compute-0; 

if (<x-3)<«N-2 && (y-2)<=(M-2)> { 
csexmin2mod3x3 * ( (x-2) %3) *3 ; 
maxdif f_compute » 

maxl3 (abs (gauss_xy_pixels [csexmin2mod3x3 + cseyminlmod3J 

- gauss_xy_middle) . maxdif f_compute) ; 

maxdif f compute = 

maxl3Tabs (gauss_xy_pixels [csexmin2mod3x3 + cseymod3J 

- gauss_xy_middle) , maxdif f_compute) ; 

maxdif f compute » 

maxl3Tabs(gaus_xy_pixelst(x%3)*3 + cseymod3] 

- gauss_xy_middle> , maxdif f_compute) ; 
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* 



cseymodB ■ -1 ; 
for ty=0; y<M+3; ++y) 

cseyminlraod3 = cseymod3; 

cseymod3 = y%3; 

cseymin2mod3 = 3-cseymod3-cseyminlmod3 ; 
for (x»0; x<N+5; ++x) { 

if (x>°4 t& x<=N+3 && y>»3 && y<=M+l) { 
csexmod3x3 = (x%3)*3; 
csexmin2raod3x3 = ((x-2)%3)*3; 
csexmin4mod3x3 * ( (x-4) %3) *3 ; 
csexmin5mod3x3 = t(x-5)%3)*3; 
if((x-5)>=l &fc (y-3)>=l) I 
if (out_compute == 255) { 

if {comp_edge_pixels [csexmin4mod3x3 +cseyrain2mod3] <comp_edge_middle) out_compute=0 

if (comp_edge_pixels [csexmin4mod3x3 +cseyminlmod3] <comp_edge_middle) out_compute*0 

if (comp_edge_pixels lcsexmin5mod3x3 +cseyminlmod3) <comp_edge_middle> out_compute=0 

if ((x-3)<=N-2 && (y-2)<=(M-2>) { 
maxdif f_compute = 

max!3 (abs(gauss_xy_pixels [csexmin2mod3x3 + cseyrainlmod3] 

- gaus9_xy_middle) , maxdif f_compute) ; 

maxdif f_compute = 

maxl3 (abs {gauss_xy_pixels [csexmin2tnod3x3 + cseymod3] 

- gauss_xy_middle) , maxdif f_compute) ; 

maxdif f_compute = 

maxl3 (abs (gaus_xy_pixels [csexmod3x3 + cseyrood3] 

- gauss_xy_middle) , maxdif f_comp\ite) ; 

)> " 
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cseymod3 * -1; 

for (y=0; y<M+3; ++y) 

cseyminlmod3 = cseymod3; 

cseymod3 = y%3; 

cseymin2mod3 = 3-cseymod3-cseyminlmod3 ; 
for (x=0; x<M+5; ++x) { 

Csexmod3x3 = (x%3)*3; 

csexmin2mod3x3 = ((x-2)*3)*3; 

csexmin4mod3x3 = ({x-4)%3)*3; 

csexminSmod3x3 = ((x-5)%3>*3; 

if (x>=4 && x<=N+3 &t y>»3 && y<=M+l) { 
if ( (x-5) >«l && (y-3) >»l) | 
' if (out compute « 255) { 

if (comp_edge_pixelslcsexmin4mod3x3 +cseymin2mod3] <comp_edge_middle) out_compute=0 

if (comp_edge_pixels [csexmin4mod3x3 +cseyminlmod3] <comp_edge_middle) out_compute=0 

if(comp_edge _pixels [csexmin5mod3x3 +cseyminlmod3] <comp_edge_middle) out_compute=0 

if ((x-3)<=N-2 && (y-2)<=(M-2)) { 
maxdif f_compute = 

maxl3 (abs (gauss_xy_pixele [csexmin2mod3x3 + cseyminltnod3 J 

- gauss_xy_middle) , maxdif f_compute) ; 

maxdif f compute = 

maxlsTabs (gauss_xy_pixels [csexmin2mod3x3 + cseymod3J 
-' gauss_xy_middle> , maxdif f_compute) ; 

maxdif f_compute = 

maxl3 (abs (gaus_xy_pixels [csexmod3x3 + cseymod3J 

- gauss_xy_middle) , maxdif f_compute) ; 

, ) '" 
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cseymod3 ■ -1 ; I 
for (y°0; y<M+3; + +y) " 
cseyminlmod3 * cseym< 




cseymod3 = y%3; 
cseymin2mod3 » 3-cseymod3-cseyminlmod3 ; 
csexmod3x3=-3 ; 
for (x=0; x<N+5; ++x) { 

csexminlmod3x3 = csexmod3x3 ,* 
csexmod3x3 = (x%3)*3; 

csexmin2mod3x3 = 9-csexmod3x3-csexminlmod3x3 ; 

if (x>=4 && x<=N+3 && y>=3 && y<=M+l> { 
if{(x-5)>=l && (y-3)>=*l> { 
if (out compute « 255) { Jja , 

if (comp_edge_pixels fcsexminlmod3x3 +cseymin2mod3] <comp_edge_middle) out_compute=0 

if (comp_edge_j>ixels [csexminlmod3x3 +cseyminlmod3] <comp_edge_middle) out_compute=0 

if (comp_edge_pixels [csexmin2mod3x3 +cseyminlmod3] <comp_edge_middle) out_compute=0 



if ((x-3)<=N-2 && (y-2)<=(M-2)) { 
maxdif f_compute = 

maxl3 tabs (gauss_xy_pixels tcsexmin2mod3x3 + cseyminlmod3] 
- gauss_xy_middle) , maxdif f_compute) ; 



maxdif f_compute = 

maxl3 (abs (gauss_xy_pixels [csexmin2mod3x3 + csey 
- gauss xy_middle), maxdif f_compute) ; 



maxdi f f _compute = 

maxl3 (abs (gaus_xy_jpixels [csexmod3x3 + cseymod3] 
- gauss xy_middle), maxdif f_compute) ; 



} 
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cseymod3 = 

for (y=0; y<M^^^+y) { 
cseyminlmod3 » cseymod3 ; 
cseymod3 = y%3; 

cseymin2mod3 = 3-cseymod3-cseyminlmod3 ; 
csexmod3x3 = - 3 ; 
cseyminlmod2 = (y-l)%2; 
cseymod2=l-cseyminlmod2 ; 
for (x=0; X<N+5; ++x) { 

csexminlmod3x3 = csexmod3x3 ; 

csexmod3x3 = {x%3)*3; 

csexmin2mod3x3 » 9-csexmod3x3-csexminlmod3x3 ; 
csexminlx2 = (x-l)*2; 
csexmin3x2 = csexminlx2-4 ; 

if (x>=3 && x<N+3 y>=2 && y<M+2) 
tmparray t (csexmin3x2+ cseymod2) %160 

+ (csexmin3x2+ cseymod2) /160*256 + 96] 
= comp_edge_pixels [csexmod3x3 

+ cseymin2mod3] = maxdif f_compute; 

if (x>= 1 && x<N+l Sl& y>=l && y<=M) 

tmparray [ (csexminlx2 + cseyminlmod2 ) %64 

+ (csexminlx2 + cseyminlmod2) /64*256] 
= gauss_xy_pixels [csexminlmod3x3 

+ cseyminlmod3] = gauss_xy_compute; 
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cseymod3 = - 1 ; 
for (y»0; y<M+3^^^^) { 

cseyminlmod3«cseymod3 ; 
cseymod3 ++; 
if(cseymod3 >= 3){cseymod3 -= 3;} 
cseymin2mod3 = 3-cseymod3 -cseyminlmod3 ; 
cseyminlmod2 « (y-l)&l; 
cseymod2=l-cseyminlmod2 
csexmod3x3= -3; 
csexx2modl 60_l_2=cseymod2 - 8 ; 
csexx2divl6O_l_2=0; 
for (x=0; x<N+5; ++x) { 

csexminlmod3x3=csexmod3x3 ; 
csexmod3 ++; 

if{csexmod3 3) {csexmod3 -= 3;} 
csexmod3x3ecsexmod3*3 ; 

csexmin2mod3x3 = 9-csexmod3x3-csexminlmod3x3; 
csexx2modl60_l_2+=2; 
csexminlx2 = (x-l)*2; 
csexmin3x2 » csexminlx2 -4 ; 

if (csexx2modl60_l_2>=160) {csexx2modl60_l_2-=160;csexx2divl60_l_2++; } 

if (x>*3 && x<N+3 && y>=2 && y<M+2) 

tmparray lcsexx2modl 6 0_1_2 + csexx2divl60_l_2*256 + 96] 
= comp_edge_pixels [csexmod3x3 

+ cseymin2mod3j - tnaxdif f_compute; 

if (x>* 1 && x<N+l && y>=l && y<^M) 

tmparray [( {csexminlx2 + cseyminlmod2) &63) 

+ ((csexminlx2 + cseyminlmod2) >>6 ) *256] 
= gauss_xy_pixels [csexminlmod3x3 

+ cseyminlmod3] = gauss_xy_compute ; 
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• 



Consider expressions in 2 different 
conditionals . 
(1) The 2 conditionals ranges are overlapping ? 





( . The expressions 

' I shows common factors ? 



No CH for these expressions 




(4)rangel + range 2 > whole range ? 




No CH for these expressions 



(5) Evaluate: 

- the ranges compared to the whole range (kl and Ic2) 

- the cost of expressions in each branch (cl. and c2) 

- and the cost after a potential CH (C_new) 

- degrees of similarity (S - (cl+c2-C_new) /C_new ) 

cl+c2 < (l+S>Mcl*kl + c2*k2) ? 



Do CH across 
conditional for 
these expressions 



No CH for these 
expressions 





Do CH across 
conditional for 
these expressions 
S=l -> <S) 
S=0 *> (3)/no 



No CH for these 
expressions 
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